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The violation of a classical Cauchy-Schwarz (CS) inequality is identified as an unequivocal sig- 
nature of spontaneous Hawking radiation in sonic black holes. This violation can be particularly 
large near the peaks in the radiation spectrum emitted from a resonant boson structure forming a 
sonic horizon. As a function of the frequency-dependent Hawking radiation intensity, we analyze 
the degree of CS violation and the maximum violation temperature for a double barrier structure 
separating two regions of subsonic and supersonic condensate flow. We also consider the case where 
the resonant sonic horizon is produced by a space-dependent contact interaction. In some cases, 
CS violation can be observed by direct atom counting in a time-of- flight experiment. We show that 
near the conventional zero-frequency peak, the decisive CS violation cannot occur. 
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The emission of Hawking radiation from the horizon of 
a black hole is an intriguing prediction of modern physics 
PQ. Due to the extremely low effective temperature, 
its detection in a cosmological context is unlikely to be 
achieved in the foreseeable future. However, it was noted 
by Unruh [2, 3 that Hawking radiation (HR) is an essen- 
tially kinematic effect that could be observed on a lab- 
oratory scale at temperatures which, while still too low, 
lie within conceivable reach. For a quantum fluid pass- 
ing through a sonic horizon (i.e., a subsonic-supersonic 
interface), it has been predicted [MTU] that, even at zero 
temperature, phonons will be emitted from the horizon 
into the subsonic region. Attempts have been made to 
observe HR in an accelerated Bose-Einstein (BE) con- 
densate [TT]. An alternative route may be provided by 
a quasi-stationary horizon, which can be achieved by al- 
lowing a confined large condensate to leak in such a way 
that the outgoing beam is dilute and fast enough to be 
supersonic [T2*l H5] . 

Hawking radiation is a fundamentally quantum phe- 
nomenon that results from the impossibility of identifying 
the vacuum of incoming quasiparticles with that of outgo- 
ing quasiparticles [14) . Specifically, the incoming vacuum 
is a squeezed state of outgoing quasiparticles. In this re- 
spect, it has been long recognized in quantum optical 
contexts [151 [Tf5] that correlation functions characteriz- 
ing the electromagnetic radiation satisfy Cauchy-Schwarz 
(CS) type inequalities that can however be violated in the 
deep quantum regime. Thus, violation of CS inequal- 
ities is generally regarded as a conclusive signature of 
quantum behavior. By contrast, detection schemes based 
on the space correlation function 6 show no signal dif- 
ference between spontaneous and thermal (stimulated) 
Hawking radiation processes [P7] . 

An additional advantage of the focus on the violation 
of CS inequalities is that it permits to distinguish the 
specific squeezed character of HR from the general prop- 
erties of coherent collective behavior. Measurements of 



space correlation functions [B], or phonon |18j or atom 
[12") intensity spectra, would not allow for such a dis- 
tinction. The reason is that the same Bogoliubov-de 
Gennes equations describe two different phenomena: lin- 
earized collective motion and quantum quasi-particle ex- 
citation. Collective motion is imprinted on the coher- 
ent (condensed) part of the wave function, which does 
not describe the squeezed zero-point dynamics of Bogoli- 
ubov quasi-particles. Importantly, we propose that spon- 
taneous HR can be tested by a CS violation involving two 
specific outgoing channels, one traveling against the flow 
in the subsonic region and the other one dragged by the 
flow on the supersonic side. CS violation due to other 
phenomena will appear in other correlation functions. 

It has recently been noted that HR could be observed 
more easily in contexts where the predicted spectrum is 
not thermal but peaked around a discrete set of frequen- 
cies [HJ [T3 • This could be the case in a sonic horizon 
formed by a double-barrier structure [T^] with an effec- 
tive potential V(x) — Z[6(x) + S(x — d)] lying between 
the subsonic and supersonic asymptotic regions. Opti- 
cal analogs can show similar peaked structures [IS]. The 
purpose of this paper is to show that, despite the re- 
maining difficulties, the violation of CS inequalities in 
Hawking radiation is comparatively much easier to ob- 
serve near the peaks characteristic of resonant radiation. 
The quest for CS violation here proposed complements 
those approaches relying on the direct detection of en- 
tanglement, as applied to inflationary cosmology [20[ . 
black holes [5T], general relativistic quantum fields |22[ . 
or black-hole (BH) analogs [23"1 124] . 

We focus on the properties of the normalized second- 
order correlation function, defined for light as [T5J [IB] 
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where is the Heisenberg operator for photon mode i, 
and the average is quantum-statistical. The correlation 
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function for classical light is obtained by removing the 
quantum average and replacing the Hcisenbcrg operators 
fit(t) by complex numbers. The following inequalities can 
be proven for classical light: 
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These inequalities are satisfied not only by any classical 
state but also by quantum thermal states at high tem- 
perature. They are also satisfied by chaotic and coherent 
states. 

The violation of any of the above inequalities is a signa- 
ture of deep quantum behavior. States violating ^ are 
said to show sub-Poissonian statistics. Violation of |3]) 
reflects anti-bunching. Expression Q is a CS inequal- 
ity. States which violate it at t — are said to exhibit 
two-mode sub-Poissonian statistics. In general, the proof 
of @ requires the system to be described by a positive 
(Glauber-Sudarshan) P function. Some quantum states 
such as two-mode squeezed states may not satisfy this 
condition and thus can violate Q. That could be the 
case in a collision between two BE condensates [25] . 

References [H [T31 El have addressed the existence 
and possible detection of HR in bosonic condensates. 
Figure [I] shows the dispersion relation of the scattering 
channels following the mode notation of Refs. [T2"l IT3"1 ITT1 
to which the reader is referred for further details. A 
particular type of Hawking radiation often considered 
in analog systems is the emission of u-out phonons as 
originated in the anomalous transmission from the <i2-in 
channel. Hereafter, the operator 7,- iQ destroys a quasi- 
particle in mode i-a, with i = u,dl,d2 and a = in, out. 
The dependence of ji a (u>) on the quasiparticle frequency 
oj will often be understood. As is conventional in fi- 
nite temperature HR setups, we assume averages are 
taken for a thermal distribution of incoming quasiparti- 
cles, so that {%i n {^)lj,m{^')) = rii(uj)Sij5(uj — uj'), where 
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FIG. 1: Dispersion relation on the subsonic (left, upstream) 
and supersonic (right, downstream) sides. The blue/red 
branches correspond to positive/negative normalization. As 
in Refs. H2lll3lfT71 d/u denotes downstream/upstream. Here, 
£ u /d denotes the asymptotic healing length, c u / d and v u / d the 
sound and flow velocities, and oj max the frequency above which 
no Hawking radiation can be generated. 



m{u)) = [exp(hQi(uj) I ksT) — and £l;(w) is the co- 
moving frequency corresponding to the mode z-in at the 
laboratory frequency uj. 

We consider the equal-time second-order correlation 
function for the outgoing quasiparticle operators 



r, 



= (7i,out7] 1 out7i,out7»,out) > 0, 



(5) 



and define io = Iy/ y/TuTjj, 6^ = - y/TuTjj, not- 
ing that the CS inequality Q is violated if and only if 

% > 1 (or Q VJ > 0) . (6) 

Thus, we may use 9ij (0jj) as a relative (absolute) figure 
of merit to quantify the degree of CS violation. 
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define the complex vector 
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element ij of the scattering matrix S characterizing the 
transition from j-in to i-out j26) and which obeying the 
pseudo- unitary condition S'r)S = rj = diag(l, 1, — 1). 
Wick's theorem allows us to write 

r d2 d2 = 2(|a d2 | 2 -I) 2 

T ud2 = \al-a d2 \ 2 + \a u \ 2 {\a d2 \ 2 -1). (7) 

Making use of ([7]), the CS inequality Q for outgoing 
quasiparticles can be rewritten as 



\al ■ a d2 \ 2 < K| 2 (|a d2 | 2 - l) , 



(8) 



which, due to the negative term within the bracket, can 
be violated some times. Interestingly, the very possibility 
of violating the CS inequality Q is a direct consequence 
of the anomalous character of the scattering process d2-in 
— > u-out, because the uj dl channel has positive/negative 
normalization. In fact, for the (normal) conversion dl o 
u, we obtain that ^ amounts to |apadi| 2 < la^pjadil 2 , 
which is always satisfied. 

The inequality ^ and pseudo- unitary lead, after a 
lengthy calculation, to the equivalent relation 

\S u d2\ 2 (l + n u + n dl + n d2 ) < 

\S d \ u \ 2 n d in d2 + \S d i d i\ 2 n u ii d2 + \S d i d2 \ 2 n u n d i 

+ \S d2d i\ 2 n u + \S d2u \ 2 n d i . (9) 

A similar inequality can be derived for the other anoma- 
lous process, dl O d2 (Andreev reflection [27]). by inter- 
changing u and dl in 

In the conventional (cj = 0) peak of the Hawking spec- 
trum \S ud2 (io)\ 2 , the scattering matrix elements diverge 
as Sij(u)) ~ with i arbitrary and j = dl,d2 . By 

contrast, Si U (uj) in the same limit (uj — > 0) tends to a 
nonzero constant. On the other hand, the only occupa- 
tion factor which diverges is n u (ui) ~ because Q u (uj) 
is the only comoving frequency that vanishes for small 
uj. From pseudo-unitarity it follows | *S' ti rf2 1 2 — |S , d2(ii| 2 = 
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IS^mI 2 — \Sd\d2\ 2 - We conclude that Q cannot be vio- 
lated in the u — ¥ region. This argument relies only on 
the presence of a constant flow connecting the subsonic 
and supersonic asymptotic regions, and not on other de- 
tails of the scattering structure. 

The inequality ([9| is however manifestly violated at 
temperature T = and w/0 provided S u d2 ^ 0, which 
further reflects the direct link between CS violation and 
Hawking radiation. Complementarily, the condition ([6]) 
is equivalent to 
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where pseudo-unitarity has again been invoked. The con 
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dition ( 10 1 is guaranteed to be satisfied because pseudo- 
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unitarity requires \Sd2d2 
nally, we note that Q d \d2 = ®ud2 
zero temperature. 

A device producing HR works like a non-degenerate 
parametric amplifier, which is known to generate squeez- 
ing from vacuum. However, sources other than vacuum 
may also generate squeezing and ultimately CS violation 
|16j . In particular, the absolute amount of CS violation 
for a given frequency often increases initially as the tem- 
perature is raised from zero, eventually reaching a maxi- 
mum and decreasing to zero at high temperatures, as can 
be guessed from Eq. ([9]). Therefore, one may wonder to 
what extent the CS violation here contemplated would 
provide conclusive evidence of spontaneous (zero-point) 
HR radiation. The answer is yes, as we argue below. 

Wick's theorem can again be invoked to write 
in terms of first-order correlation functions such as 
(7i,out7j,out) or (7i,out7j,out) ; here generically referred to 
as pij . The thermal contribution is pff = p^ — pij , with 
p®j the zero-temperature value (average over incoming 
vacuum) . If one neglects the zero-point contributions by 
approximating pij ~ pfj, then one arrives at a modified 
version of ([9| where only the terms quadratic in the rij 's 
survive. The resulting inequality is always satisfied. We 
conclude that CS violation requires vacuum fluctuations. 

We may define T v as the highest temperature at which 
([9]) is violated is satisfied] and try to identify some 
trends in its behavior. We note that, if the upstream 
current is small (kinetic contribution to total chemical 
potential on subsonic side is small), then in the central 
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regions of the Figs. 



(w ~ w max /2) the comoving 
frequencies of the various channels are comparable (call 
them CI') and all smaller than (typically a fraction of) 
the comoving chemical potential on the subsonic side, 
CI' < p. The resonant BH configurations, first discussed 
in [i2j , have the important advantage of displaying HR 
peaks at nonzero frequencies which generally lie in the 
interesting region near cj max /2. 

If the conversion from negative to positive normal- 
ization is weak, then | <Srf2d2 1 is close to unity and from 



pseudo-unitarity it can be proven that all the S-matrix 
elements appearing in ^ are small except for the rela- 
tion \S dlu \ 2 + \S dld i\ 2 ^ 1. If L = -log|S„ d2 | > 1 we 
obtain a low violation temperature, T v ~ CI' /L -C p- 

Conversely, we find that for the highest HR peaks 
(where |S ud2 |»l), then either |S„ d2 | ^ \S d 2u\ > \S d 2di\ 
or \S ud2 \ ^ \S d 2di\ > \S d 2u\- At the same time, \S ud2 \ > 
I jSdij I for all j, which implies that the violation tempera- 
ture T„ ~ Cl'\S ud 2\ 2 /Sj 1 > p, where S d i = maxj{\Sdij\}- 
Thus, T v {uj) is expected to be large near the peak fre- 
quency wo- Non-resonant structures can also show CS 
violation at nonzero temperatures. However, even if the 
relative violation (as measured by u d2) is significant, the 
absolute amount of violation (as measured by O u d 2 ) is 
negligible, as discussed later. 

Another important advantage of resonant peaks is 
that, at their relatively high frequency, the phononic 
signal is approximately proportional to the atomic sig- 
nal, which can be directly measured in a time-of-flight 
(TOF) experiment. Assuming that one can separate the 
subsonic and supersonic TOF signals, near the peak at 
lj = ujq, and neglecting finite-size effects, the atom oper- 
ators can be approximated as: 

(11) 
(12) 
(13) 

where c u /d{k) annihilates an atom of momentum k on 
side u/d. Here, fc, )0 ut is the comoving quasiparticle mo- 
mentum at the laboratory frequency u in the i-out chan- 
nel, and q u /d is the condensate momentum per atom on 
each side. The proportionality factors not shown in (11- 
13 1 cancel in the CS violation condition (|6| for 8 ud2 and 

#(21(22 • 

For the approximation ( 12p3 1 to apply, it is necessary 
that, for scattering processes with output on the super- 
sonic side, the CS- violating S'-matrix elements must dom- 
inate over the S'-matrix elements at different frequencies 
but with the same laboratory momentum pi 7 implicitly 
defined in ( TT|p~3 ). On the subsonic side, the approxima- 
tion (111 requires the peak to be large, \S u d2(^o)\ 2 ^ 1) 
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because it has to stand out above the background of the 
depletion cloud [T2"] . 

From (10 1 and the pseudo-unitarity relation | *S r u .d2 1 2 + 
|Sdi(22| 2 + 1 = \Sd2d2\ 2 , it may appear that a shortcoming 
of a large peak is its small relative degree of CS-violation 
(as measured by 9 u d2)- However, this does not imply 
that the experimental signal is necessarily small. Quite 
the opposite, the absolute amount of violation (as mea- 
sured by O u( 2 2 ) can be quite large, as (10 1 directly re- 
veals. A similar analysis can be performed for the other 
anomalous process, dl d2. 

Whether or not the conditions for the approximation 
( 11|13 | are satisfied, we can define, in analogy to ([5]), the 
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FIG. 2: Hawking radiation for a condensate leaking through 
a double barrier structure like that studied in Ref. I12I The 
strength of both delta barriers is Z = 2.2h 2 /m£ u . They 
are separated by a distance d — 3.62£ u . The flow is such 
that q u £,u = 0.01 and o; max = 0.99/i/ft. Solid blue: zero- 
temperature Hawking radiation spectrum \S u d2(u)\ 2 ■ Dashed 
red: 9 u d2(w) at temperature T = /u/fes, where /x is the comov- 
ing chemical potential on the subsonic side. Dotted green: 
maximum violation temperature T v (oj) in units of /i/fce. Not 
shown in the figure, T v (uj) rises up to T v (wo) — 21/i/fcs, 
where |Sw2(wo)| 2 — 8. Left inset: Zoom of the peak region, 
with T v (u) removed and z u( n(w) (relative atom CS violation) 
added (dotted purple). Right inset: same as left inset; it 
shows Q u di(y>) (dashed-dotted brown) and Z u di(ui) (dotted 
purple), which measure the absolute amount of CS violation 
in the phonon and atom signals. 



atomic correlation functions: 

G„ d2 (w) = (ci(p u )c d (p d2 )c d (p d2 )c u (p u )), (14) 
and similarly for the other Gy, where i,j — u,dl,d2, 



and Pi(u>) is defined in ( TTfl3 1. It can be shown that a 
sufficient condition for ([6j) is: 



> 1 



Z tj > . 



(15) 



where Zij = (/,, \ (i, ,(•', , and Zij = Gij — \J Gu Gjj . 

In Fig. [2] we plot, for a double delta-barrier structure, 
the zero temperature HR spectrum \S ud 2{io)\ 2 , together 
with 6 u d%{uS) (at ksT = fi) and the violation tempera- 
ture T v (iS). The left inset magnifies the peak region and 
includes z ud2 {iS)^ which measures the relative amount of 
CS violation in the atomic signal. The right inset shows 
©ud2(w) and Z ud2 (oj), i.e., the absolute amount of CS 
violation by the phonon and atom signals. The latter 
can be directly measured in a TOF experiment. These 
graphs indicate that the considered structure is a promis- 
ing scenario for the unambiguous detection of sponta- 
neous Hawking radiation. 

The inset of Fig. [3] shows the corresponding curves for 
a single delta-barrier black-hole configuration. The HR 
spectrum is unstructured, with a single peak at co=0. 
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FIG. 3: Same as Fig. [2] for a structure without barriers but 
with two sharp variations in the local speed of sound, which 
takes the successive values c„, 0.5c u , 0.2c„. The intermediate 
region has a length 12£„. The flow is such that q u £u = 0.7 
and ui max = 0.21/x/ft. This is a resonant generalization of the 
model first studied in Ref. [6| where only one sound-speed 
discontinuity was considered. Inset: Same curves for a setup 
with a single delta barrier of strength Z = 0.62ft 2 /m£„ and a 
uniform interaction, with flow q u £,u = 0.3 and w ma x = 0.6/i/fi. 
Atomic CS violation is not found for these structures. 



The CS inequality can be violated at the relatively high 
temperature T — fi/ks within a considerable frequency 
range. However, due to the smallness of IS^^I, the 
atomic signal for a given momentum includes contribu- 
tions from different phonon frequencies and from the de- 



pletion cloud, which implies that the approximation (11 
13 1 cannot be used. Moreover, the smallness of the abso- 



lute phonon violation signal, together with the absence 
of atomic CS violation, suggests that non-resonant struc- 
tures are not good candidates for the observation of CS 
violation. 

The main Fig. [3] shows the same curves as in Fig. 
[2] but for a different type of resonant structure where 
spatial dependence of the interatomic contact interaction 
is tuned to produce step-like variations in the local speed 
of sound. This model is inspired in the setup proposed in 
Ref. O where the subsonic-supersonic interface is induced 
by a single step in the profile of the speed of sound. The 
absolute CS violation is considerably smaller here than 
in Fig. [51 but could still be observable. The sufficient 



atom inequality ( 15 ) is not satisfied 



In conclusion, we find that CS violation in the strongly 
peaked Hawking radiation emitted by a double-barrier 
sonic black-hole may provide a convenient route to the 
unambiguous observation of the zero-point contribution. 
In some setups, the violation of CS inequalities is large 
enough to be detectable by the direct observation of 
second-order correlation functions in an atom time-of- 
flight of experiment. These advantages are absent in 
non-resonant structures such as that formed by a single- 
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barrier sonic black-hole. In particular, the relevant 
CS violation cannot occur near the conventional, zero- 
frequency Hawking radiation peak universally shown by 
all one-dimensional BH structures. 
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